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A description of phonon dynamics and charge response of the n-doped high-temperature super- 
conductor (HTSC) NdCuO is presented based upon a microscopic modehng of the electronic den- 
sity response. This is accomplished starting from the insulating state via the underdoped strange 
metallic to the more conventional metallic state by consecutive orbital selective incompressibility- 
compressibility transitions in terms of strict sum rules for the charge response. The approach 
proposed in this work for the n-doped HTSC's modifies the modeling recently applied to the p- 
doped compounds and expresses an electron-hole asymmetry introduced by doping. A qualitative 
physical picture consistent with our modeling of the electronic state in the cuprates is given in 
which a sufficiently broad set of orbital degrees of freedom, i.e. Cu3d/4s and 02p is essential. 
Within the framework of linear response theory we calculate full phonon dispersion curves in the 
different phases. In particular, the strongly doping dependent anomalous high-frequency oxygen 
bond-stretching modes (OBSM) found experimentally in the p-doped HTSC's and also recently for 
n-doped Ndi.85Ceo.i5Cu04 are investigated and compared with experimental results from inelastic 
neutron scattering and inelastic x-ray scattering, respectively. We calculate an anticrossing scenario 
for the OBSM in n-doped NdCuO which is absent in the case of p-doped LaCuO and relate it to 
the different crystal structure. Phonon-induced electronic charge redistributions of the anomalous 
OBSM due to nonlocal electron-phonon interaction effects of charge-fluctuation type giving reason 
to dynamic stripes are also studied. Finally, calculations of a characteristic rearrangement of the 
phonon density of states across the insulator-metal transition are presented and a comparison with 
experimental results has been accomplished. 
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I. INTRODUCTION 

There is increasing evidence that the electron-phonon 
coupling is strong in the cuprate based high-temperature 
superconductors (HTSC's) and phonons might play an 
important role for the electron dynamics in the HTSC's, 
see e.g. Refs. miiiiilli. It has been known 
from the experiments for some time that the frequen- 
cies of the high-frequency oxygen bond-stretching modes 
(OBSM) are strongly renormalized (softened) upon dop- 
irig in t he p- dope d cuprates investigated so far, Refs. 
illi^lilllililHlil. More recently, an anomalous soft- 
ening behaviour of the OBSM, not present in the un- 
doped insulating state, also has been observed in the n- 
doped metallic state of NdCuO, Refs. These 
findings support the generic nature of the anomalies in 
all HTSC's. The frequency renormalization and a cor- 
responding increase of the phonon linewidths point to a 
strong coupling of these phonons to the charge carriers. 
Theoretical calculations for the p-doped cuprates have 
shown that the strong nonlocal electron-phonon interac- 
tion (EPI) leading to the softening of the OBSM is due 
to a specific screening effect in terms of a local polariza- 
tion created by metallic charge fluctuations (CP's) on the 
outer shells of the Cu an Oxy ions in the CuO-plane, see 
e.g. Refs. [HllOIlll. 

This type of screening in a strongly inhomogeneous 
electronic system induced by a density distribution which 



mainly is concentrated at the ions in the CuO-plane leads 
for the OBSM in the p-doped cuprate a^^i^^ and also for 
n-doped NdCuO, as shown in this work, to an anoma- 
lous behaviour of the dispersion. The specific screening 
yields a downward dispersion of the longitudinal branch 
along the (1,0,0) and the (1,1,0) direction falling in the 
(1,0,0) direction in LaCuO and partly also in NdCuO 
even below the transverse branch. Such a behaviour 
cannot be understood within a typical lattice dynamical 
model, like e.g. the shell model, including a homoge- 
neous electrongas screening that is diagonal in reciprocal 
space. The latter type of dielectric response always gives 
an increasing dispersion when passing from the center of 
the Brillouin zone (BZ) into the zone, because of incom- 
plete screening of the changes of the Coulomb potential 
at shorter distance. Such a behaviour is found in our 
calculations for the dispersion of the branches for modes 
with polar character at lower energy. 

Contrary, the anomalous softening of the OBSM lead- 
ing to a decreasing dispersion is a result of an "over- 
screening" at shorter wave lengths of the changes of the 
Coulomb potential generated by the motion of the ions in 
these particular modes. Physically, this "over-screening" 
in the OBSM giving rise to bond length modulations has 
been shown by our calculations to be due to nonlocally 
excited ionic CP's localized at the Cu and Oxy sites lead- 
ing to a (dynamic) charge ordering in form of localized 
stripes of alternating sign in the CuO plane. Expressed in 
terms of the density response matrix in reciprocal space 
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such strong nonlocal EPI effects and "over-screening" 
of the OBSM is related to their off- diagonal elements 
(local-field effect) being apparently very important in the 
cuprates as revealed by the anomalous softening. 

It is interesting to note in this context, that an "over- 
screening" effect is not found in the in-plane OBSM of 
metallic SrRuO, a low-tcmpcraturc. Tc=1.5K, (presum- 
ably) triplet superconductor, structurally isomorphic to 
LaCuO, see Ref. H^- This points to a more homogeneous 
electrongas-like screening mechanism in this material and 
correspondingly to a reduction of the strong nonlocal EPI 
found in the cuprates. The different screening in SrRuO 
as compared with the cuprates should be associated with 
a more extended Ru4(i state as compared to the localized, 
stronger correlated Cu3c? state and a strong hybridiza- 
tion between Ru4d and 02p. Furthermore, SrRuO is not 
a low-carrier metal and has a large density of states at 
the Fermi energy. 

In our microscopic theoretical description of the charge 
response, the EPI and lattice dynamics of the HTSC's 
in the framework of the linear response approach wc 
focus on the specific features of their solid-state chem- 
istry which is necessary for a realistic description of the 
cuprates since the phenomena observed in the HTSC's 
hardly occur in the other known compounds. 

The strong ionic nature of the cuprates is modeled in 
our approach by an ab initio rigid ion model (RIM) , used 
as a reference frame for the local ionic, rigid charge re- 
sponse and the EPI, respectively. The important and 
characteristic nonlocal, nonrigid contribution to the elec- 
tronic density response and the EPI in the HTSC's is cal- 
culated in terms of microscopically well defined CF's and 
dipolc fluctuations (DF's). The CF's and DF's arc ex- 
citable on the electronic shells of the ions in the crystal. 
Additionally, covalent metallic features of bonding are 
approximatively taken into accoun t^^'^'^ . In the calcula- 
tions a sufficiently broad set of orbital degrees of freedom 
(CuSd, Cu4s and 02p), the full three-dimensional long- 
ranged Coulomb interaction as well as short-ranged local 
repulsions of the electrons, in particular the important 
on-site repulsion mediated by the localized CuSd orbitals, 
is considered quantitatively. All the resulting couplings 
arising in the dynamical matrix and the EPI are micro- 
scopically well defined and can be calculated. The results 
obtained for the p-doped materials so far agree well with 
the experimental phonon dispersion. For a recent review, 
see Ref. [11. 

Our model for the (irreducible) electronic polarizabil- 
ity matrix, H^k' , specifying the kinetic part of the charge 
rcsponsei2i2i, is defined in a localized basis by the orbital 
degrees of freedom k, k' in the elementary cell of the crys- 
tal. In the calculations for the HTSC's Cu3d, Cu4s and 
02p orbitals in the CuO plane are taken into account. 

In the absence of a rigorous quantitative description 
of the electronic state of the HTSC's our model for H^k' 
has been designed for the different phases to be consis- 
tent with the corresponding charge response in the long- 
wavelength limit. This is achieved by fulfilling rigorous 



sum rules for the electronic density response in terms of 
the polarizability matri5*22i^. These sum rules, see also 
Sec. ini can be considered as orbital resolved closed forms 
to represent the change of the charge response across a 
metal-insulator transition in terms of the electronic po- 
larizability or the compressibility of the electronic sys- 
tem, respectively, as a primary tool to characterize the 
corresponding ground state. 

Hitherto such a modeling of the electronic state of the 
cuprates only has been applied to the p-doped materials, 
see Refs. [l9, 21,23 for a detailed discussion. We start 
from the insulating state of the HTSC's where both the 
Cu3d and 02p orbitals are taken to be incompressible, 
consistent with a gap in the charge excitation spectrum, 
and the delocalized Cu4s orbitals are neglected. 

For the underdoped " strange metallic" phase of the p- 
doped HTSC's, where the description of the low-energy 
excitations remains a theoretical challenge, we have pro- 
posed a pseudogap model for the electronic density re- 
sponse in the pasii^^. It describes the localization- 
delocalization transition in the electronic structure in 
terms of a qualitatively different compressibility of the 
Cu3d and 02p states, respectively, and admits that the 
single-particle density of states is suppressed at the Fermi 
energy for the strongly localized and correlated Cu3d 
states because of the cost in energy from hopping of the 
charge carriers to the Cu-sites. This is achieved by in- 
stalling an insulator-like, incompressible charge response 
for the localized Cu3d states via the sum rules, but not 
so for the more delocalized 02p orbitals where the holes 
predominantly are injected to in p-type cuprates. For the 
02p orbitals a compressible, metallic charge response is 
allowed with a renormalized partial density of states at 
the Fermi energy sp consistent with a gain in kinetic 
energy by delocalization. 

The blocking of the metallic charge response at the 
Cu sites in the underdoped phase is lifted in our mod- 
eling of the optimally and overdoped metallic phase by 
allowing all orbitals, i.e. Cu3d, Cu4s and 02p, to be- 
come compressible, metallic with a nonvanishing partial 
density of states at ep for all orbitals considered. In a fc- 
space picture a large Fermi surface (FS) as found, e.g., in 
LDA-calculations, can now be expected to develop from 
a gapped FS of the underdoped phase seen in the exper- 
iments by angle- resolved photoemsission (ARPES). 

Thus, according to the modeling, wc have a crossover 
in the electronic properties corresponding to a qualitative 
change of the ground state to a Fermi-liquid like state in 
the overdoped material. Apparently, the behaviour is 
quite different on the two sides of the crossover point, 
and there is no phase transition with symmetry breaking 
assumed in our modeling. 

For the n-doped materials we investigate quantita- 
tively in this work a model proposed recently in Ref. 
[T9I that modifies the sequence of the orbital selective 
incompressibility-compressibility transitions. In this way 
an electron-hole asymmetry is introduced in our model- 
ing of the insulator-metal transition, simulating the dif- 
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fercnt character of hole and electron doping. While in the 
p-doped HTSC's hole carriers are brought to the 02p or- 
bitals making the latter compressible in our picture, in 
the n-doped materials, like Nd2-i^Cea;Cu04, electron car- 
riers are brought into orbitals with Cu3(i-4s character at 
the Cu sites. Thus, we model with the help of the sum 
rules an underdoped n-type material in our approach in 
terms of a metallic, compressible charge response at the 
Cu sites (Cu3(i) with an increasing Cu4s component upon 
doping and a localized incompressible, insulator-like re- 
sponse of the 02p orbitals at the O^y sublattices of the 
CuO plane. 

Such a modeling also may be consistent with a sta- 
bilization of antiferromagnetic spin correlations, because 
an insulator-like localization of the incompressible 02p 
orbitals should be favourable for superexchange. Experi- 
mentally it is found that the phase diagram is asymmetric 
with respect to electron and hole doping, and for the n- 
type materials the antiferromagnetic phase extends much 
further with doping. From our results for the electronic 
CF polarizability of the oxygen one might speculate that 
this quantity should play an indirect role for antiferro- 
magnetism and superconductivity in the cuprates. More- 
over, the incompressible response for 02p also generates a 
pseudogap phenomenon attenuating the density of states 

Then, in the optimally doped state a crossover to a 
metallic charge response at the formerly incompressible 
02p orbitals is admitted and so all orbitals, i.e. Cu3d, 
Cu4s and 02p are modeled as compressible by allowing 
the 02p states to become compressible too, like in the p- 
doped case. Correspondingly a change of the FS topology 
can be expected and Fermi-liquid behaviour should revive 
in the system. 

Note in this context that from ARPES experi- 
ments, Refs. [25II26I of the n-type superconductor 
Nd2-2^Cea;Cu04 it is found, that at low doping the FS 
is a small electron pocket (with volume ~ x) centered 
at (tt, 0, 0). Further doping finally leads to a large LDA- 
like FS of the hole-type (volume ~ 1 -I- x) centered at 
(tt, 7r,0). A change of FS topology upon doping, differ- 
ent from the n-doped material, also has been observed 
in the p-doped HTSC's. Here, at low doping ARPES 
experiments reveal so-called " Fermi surface arcs" (which 
also could be elongated Fermi surface pockets at nodal 
positions). With further doping in the optimally doped 
and overdoped state a large LDA-likc FS of the hole-type 
forms as in the n-doped case, Refs. [27ll28l . 

The article is organized as follows. In Sec. [TTl ele- 
ments of the theory and modeling are reviewed to set 
the frame. Section [HT] presents our calculated results of 
the phonon dispersion of NdCuO across the insulator- 
metal transition taking CF's and DF's in the electronic 
density response into account. In Sec. IIVI the possibly 
generic phonon anomalies and the anticrossing behaviour 
is discussed, while Sec. |V] deals with the phonon-induced 
charge response of the OBSM. In Sec. IVII the phonon 
density of states for insulating and metallic NdCuO is 



calculated and a comparison of the redistribution of spec- 
tral weight is provided. Finally, a summary of the paper 
is presented in Sec. IVIII and the conclusions are given. 



II. THEORY AND MODELING 

From a general point of view our treatment of the elec- 
tronic density response and lattice dynamics in terms of 
DF's and CF's can be considered as a microscopic imple- 
mentation of the phenomenological dipole-shell model or 
the charge-fluctuation models, respectively. For a gen- 
eral formulation of phenomenological models for lattice 
dynamics that use localized electronic variables as adia- 
batic degrees of freedom, see for example Ref. [2^ This 
formulation covers shell models, bond-charge models and 
charge-fluctuation models. While in this approach the 
coupling coefhcients are treated as empirical fitting pa- 
rameters, the essential point in our microscopic scheme 
is that all the couplings can be calculated. 

In the following a survey of the theory and modeling 
is presented. A detailed description can be found in Ref. 
I20I and in particular in Ref. |30| where the calculation of 
the coupling parameters is presented. 

The local part of the electronic charge response and 
the EPI is approximated in the spirit of the Quasi-Ion 
concept^ by an ab initio rigid ion model (RIM) taking 
into account covalent ion softening in terms of (static) ef- 
fective ionic charges calculated from a tight-binding anal- 
ysis (TEA). In addition, scaling of the short-ranged part 
of certain pair potentials between the ions is performed 
to simulate further covalence effects in the calculation 
in such a way that the energy-minimized structure is as 
close as possible to the experimental one^. Structure 
optimization and energy minimization is very important 
for a reliable calculation of the phonon dynamics through 
the dynamical matrix. 

The RIM with the corrections just mentioned then 
serves as an unbiased reference system for the descrip- 
tion of the HTSC's and can be considered as a first ap- 
proximation for the insulating state of these compounds. 
Starting with such an unprejudiced rigid reference sys- 
tem nonlocal, nonrigid electronic polarization processes 
are introduced in form of more or less localized electronic 
charge-fluctuations (CF's) at the outer shells of the ions. 
Especially in the metallic state of the HTSC's the lat- 
ter dominate the nonlocal contribution of the electronic 
density response and the EPI and are particularly impor- 
tant in the CuO planes. In addition, anisotropic dipole- 
fluctuations (DF's) are admitted in our approach, Refs. 
[23l30l . which prove to be specificly of importance for the 
ions in the ionic layers mediating the dielectric coupling 
and for the polar modes. Thus, the basic variable of our 
model is the ionic density which is given in the perturbed 
state by 

Pair, Qx, Pa) = P^!:(r) +^ QaP^M +Pa • rPair)- (1) 

X 
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is the density of the unperturbed ion, as used in the 
RIM, locahzed at the sublattice a of the crystal and mov- 
ing rigidly with the latter under displacement. The Q\ 
and describe the amplitudes and the form factors 
of the CF's and the last term in Eq. ([T|) represents the 
dipolar deformation of an ion a with amplitude (dipole 
moment) p„ and a radial density distribution pJJ as form 
factor, r denotes the unit vector in the direction of 
r. The p'^ are approximated by a spherical average 
of the orbital densities of the ionic shells calculated in 
LDA taking self-interaction effects (SIC) into account. 
The dipole density pj^ is obtained from a modified Stern- 
heimer method in the framework of LDA-SIC, Ref. [10. 
All SIC-calculations are performed for the average spheri- 
cal shell in the orbital-averaged form according to Perdew 
and Zungcr, Ref. [H. For the correlation part of the en- 
ergy per electron e the paramctrization given in Ref. [s^ 
has been used. 

The total energy of the crystal is obtained by assum- 
ing that the density can be approximated by a superpo- 
sition of overlapping densities p^,. The p^ in Eq. (U]) are 
also calculated within LDA-SIC taking environments ef- 
fects, via a Watson sphere potential and the calculated 
static effective charges of the ions into account. The 
Watson sphere method is only used for the oxygen ions 
and the depth of the Watson sphere potential is set as 
the Madelung potential at the corresponding site. Such 
an approximation holds well in the HTSC's^i^i. As a 
general rule, partial covalence reduces the amplitude of 
the static effective charges in mixed ionic-covalcnt com- 
pounds such as the HTSC's, because the charge transfer 
from the cations to the anions is not complete as in the 
entirely ionic case. Finally, applying the pair-potential 
approximation we get for the total energy 



i?(i?,C)=E^a(0 



(a,Q)5>^(b,/3) 



(2) 

The energy E depends on both the configuration of the 
ions {i?} and the electronic (charge) degrees of freedom 
(EDF) {(^} of the charge density, i.e., {Qa} and {Pa} in 
Eq. ([1]). E'^ arc the energies of the single ions, a, b 
denote the elementary cells and a, /3 the corresponding 
sublattices. The second term in Eq. ^ is the interaction 
energy of the system, expressed in terms of anisotropic 
pair interactions ^ap. Both E^ and hi general de- 
pend upon C via pa in Eq. ([T]). The pair potentials in 
Eq. (12) can be separated into long-ranged Coulomb con- 
tributions and short-ranged terms as follows: 

$„/3(R,C) =— ^ {2aPp + 2/3Pa) • + ^3 

-3 (P"-^^f ■P^U $.,(R,C), (3) 



$„/5(R, C) ^K^U^iR, C) + A>[/„(R, C) 
+ W„^(R,C)+G<,0(R,C)- 



(4) 



The first term in Eq. ([3]) describes the long-ranged ion- 
ion, the second the dipole-ion and the third and fourth 
term the dipole-dipole interaction. Za and Zp are the 
variable charges of the ions in case the CF's are excited. 
The latter reduce to the ionic charges for rigid ions. Ka 
and Kp are the charges of the ion cores. The remaining 
term in Eq. ([3|) given in Eq. (jj]) represents the short- 
ranged interactions. Detailed expressions for these in- 
teractions and the procedure for the calculation can be 
found in Ref. Isol . 

From the adiabatic condition 



= 



(5) 



an expression for the atomic force constants, and accord- 
ingly the dynamical matrix in harmonic approximation 
can be derived: 



RIM 



^[Br(q)]*[C-l(q)]^^,B;''(q) 



(6) 



The first term on the right hand side denotes the contri- 
bution from the RIM. Ma, Mp are the masses of the ions 
and q is a wave vector from the first BZ. 

The quantities B(q) and C(q) in Eq. ([6|) represent the 
Fourier transforms of the electronic coupling coefficients 
as calculated from the energy in Eq. ([2]). or the pair 
potentials in Eqs. ([3]), (H]), respectively: 



Bab 

K/3 



ab 



(7) 
(8) 



K denotes the EDF (CF and DF in the present model, 
see Eq. ([1])) in an elementary cell. The B coefficients 
describe the coupling between the EDF and the dis- 
placed ions (bare electron- phonon coupling), and the co- 
efficients C determine the interaction between the EDF. 
The charge and dipole fluctuations Qx and p^ do not 
appear explicitly in the dynamical matrix. Only the 
derivatives of the energy with respect to the EDF are 
needed in order to calculate the coefficients B and C, 
respectively, and these derivatives have to be performed 
at Q\ = 0, Pa = 0. The phonon frequencies uJaiq) and 
the corresponding eigenvectors e"(qcr) of the modes (qcr) 
are obtained from the secular equation for the dynamical 
matrix in Eq. ([6]), i.e. 



5]t:f(q)ef(q)=u.2(q)er(q). 



(9) 



Equations ([6])-(l9]) are generally valid and, in particular, 
are independent of the specific model for the decompo- 
sition of the perturbed density in Eq. ^ and the pair 



5 



approximation Eq. ([2]) for the energy. The lenghty de- 
tails of the calculation of the coupling coefficients B and 
C cannot be reviewed in this paper. They are presented 
in Ref. [sol . In this context we remark that the coupling 
matrix C^^' (q) of the EDF-EDF interaction, whose in- 
verse appears in Eq. ([6]) for the dynamical matrix, can 
be written in matrix notation as 



V. 



(10) 



contains the kinetic part to the interaction C and V 
the Hartree and exchange-correlation contribution. C^^ 
needed for the dynamical matrix and the EPI is closely 
related to the (linear) density response function (ma- 
trix) and to the inverse dielectric function (matrix) e~^, 
respectively. Only very few attempts have been made 
to calculate the phonon dispersion and the EPI of the 
HTSC's using the linear response method in form of den- 
sity functional perturbation theory (DFPT) within LDA, 
Refs. [mill. These calculations correspond to calcu- 
lating n and V in DFT-LDA and for the metallic state 
only. On the other hand, in our microscopic modeling 
DFT-LDA-SIC calculations arc performed for the vari- 
ous densities in Eq. ([1]) in order to obtain the coupling 
coefficients B and V . Including SIC is particularly im- 
portant for localized orbitals such as CuSd in the HTSC's. 
SIC as a correction for a single particle term is not a cor- 
relation effect, which per definition cannot be described 
in a single particle theory, but SIC is important for con- 
tracting in particular the localized CuSd orbitals. Our 
theoretical results for the phonon dispersion of p-doped 
materials, Refs. [T9ll2Tll23 . which compare well with the 
experiments, demonstrate that the approximative calcu- 
lation of the coupling coefhcients in our approach is suf- 
ficient, even for the localized CuSd states. Written in 
matrix notation we get for the density response matrix 
the relation 



c- 



n(i + Fn)-^ = ne-\ 



1 + yn. (11) 



The CF-CF submatrix of the matrix 11 can approxima- 
tively be calculated from a TBA of a single particle elec- 
tronic bandstructurc. In this case the static electronic 
polarizability 11 in tight-binding representation reads 



n,«,(q,^ = o) = -^ 



/„.(k + q)-/„(k) 
N ^^i;„.(k + q)-^„(k) 



X [c:„(k)c™,(k + q)][<,„(k)c,,„,(k + q)]*. (12) 



/, E and c in Eq. (|12p are the occupation numbers, the 
single-particle energies and the expansion coefficients of 
the Bloch functions in terms of tight-binding functions. 
The self-consistent change of an EDF on an ion induced 
by a phonon mode (qc) with frequency 0Ja{(\) and eigen- 
vector e"(qCT) can be derived in the form 



'5C(q^) 



^X'^"(q)u„(qa) 
5a(qa)e*'^^^ 



(13) 



with the displacement of the ions 



<(qa) 



2i\/„cj^(q) 



1/2 

e"(qa)e*'i^'' = u„(qcr)e*i^°. 

(14) 

The self-consistent response per unit displacement of the 
EDF in Eq. (fT51) is calculated in linear response theory 
as 



X(q) = n(q)£-l(q)B(q) = C-l(q)B(q). 



(15) 



A measure of the strength of the EPI for a certain 
phonon mode (qcr) is provided by the change of the self- 
consistent potential in the crystal felt by an electron at 
some space point r in this mode, i.e. (514ff(r, qcr). Aver- 
aging this quantity with the corresponding density form 
factor /^^(r — RJ^) at the EDF located at R)^, we obtain 



(16) 



This gives a measure for the strength of the EPI in the 
mode (qcr) mediated by the EDF considered. For an ex- 
pression of (5V)f (qcr) in terms of the coupling coefficients 
given in Eqs. ([7]) and ([5]), see Ref. [l^. From our cal- 
culations for LaCuO large values for (5V^(qcr) are found, 
in particular for the OBSM phonon anomalies and even 
larger for the nonadiabatic c-axis phonons in the metallic 
phase, mixing with a plasmon. 

A nonadiabatic approach is necessary for a description 
of the interlayer phonons and the charge response within 
a small region around the c axis, Refs. [Tgllssl . 

Finally, it may be useful to outline, how the perturbed 
charge density in Eq. ([1]) is related to the Fourier trans- 
form of the response quantity X(q) in Eq. (|15p and thus 
to n and s describing the dielectric response of a specific 
material. 

The change of the electronic density p{r), depending 
on the ionic configuration {R} and the EDF {C}, induced 
by a unit displacement of an ion at R^ is given by the 
vector field 



PS(r) 



dp(r) _ dp{r) 



OR- 



V dpjr) dQ 



(17) 



By differentiating the adiabatic condition in Eq. ([5|) with 
respect to R|^ and using the Hohenberg-Kohn functional 
of DFT for the energy the second term in Eq. (|17|1 can 
be put into the ioriv^ 



(18) 



b K 



where the form factors of the EDF's k, arc defined by 



_ dp{r) 



(19) 



and X° ^ is the Fourier transform of X*^" (q) in Eq. ([Te 
Equation ^T7\ gives a decomposition of the charge re- 
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sponsc into an explicit (rigid) change of the density, de- 
fined by 



PS.r(r) ^ 



dp{v) 



(20) 



and the distortion contribution Pi^j(r) from Eq. (fTS)) . 
Physically, the explicit contribution Pj^ ^ (r) describes the 
part of the density which rigidly follows the motion of the 
ions and represents the local EPI effects. On the other 
hand, ^{r) describes the way, how the charge density 
is distorted in response to the displacement of an ion. It 
represents the nonlocal EPI effects. 

Such an unique decomposition of the perturbed den- 
sity into a rigid and nonrigid (distortion) contribution 
has been shown to be rigorously valid in linear response 
theory, see the discussion of the Quasi-Ion concept in Ref. 
I31I . In case of a diagonal response function (matrix) in re- 
ciprocal space the distortion contribution to the charge 
response vanishes and only the rigid Quasi-Ion (Pseu- 
doatom) contribution related to local EPI survives^. Es- 
pecially, this is true for homogeneous electrongas screen- 
ing. On the other hand, in a strongly inhomogeneous 
electronic system like the cuprates the density matrix in 
reciprocal space, is of course, nondiagonal and appar- 
ently there are important distortion contributions asso- 
ciated with the strong nonlocal EPI effects which add to 
the charge response. In our modeling these distortion 
effects are described in terms of ionic CF's and DF's, re- 
spectively. This is an adequate approximation because of 
the specific solid-state chemistry of these compounds and 
has been confirmed so far by our calculations for p-doped 
materials and also for n-doped cuprates in this work. 

In the present work the rigid part is approximated by 
the rigid- ion density p° in Eq. ([1]), which works well for 
materials with a strong component of ionic binding and 
the distortion part by the two remaining terms with the 
corresponding form factors of the CF's and DF's as given 
in this equation. The latter yield the shape of the change 
in the density related to the CF's and DF's, respectively. 

A central problem of superconductivity in the cuprates 
is as to how the metallic state evolves from the Mott- 
insulator upon doping. In particular, there is no con- 
sense for a realistic description of the underdoped region. 
Thus, there remains the important question as to how to 
discriminate between the charge response of the metallic 
and insulating state of the HTSC's, respectively. Specif- 
ically, the latter cannot be obtained for example within 
the LDA and a realistic quantitative description to cal- 
culate the irreducible (proper) polarization part nKK'(q) 
for the HTSC's is not available. However, a general crite- 
rion which necessarily requires a multi-orbital approach 
follows from the different analytic behaviour of the or- 
bital resolved polarizability in the longwavelength limit 
(q 0) in both phase a^^i^" . 

In the metallic phase the electronic partial density of 
states (PDOS) at the Fermi level Z^{eF) is related to the 



polarizability matrix for (q 0) according to 
^H,,,(q^O) = Z,(£F), 



(21) 



and the total density of states (DOS) at energy e is given 

by 



Z(£)=^Z„(e). 



(22) 



On the other hand, for the insulating state we obtain the 
sum rules 



and 



5^n,„.(q^0) = O(g) 



^H«,,(q->0) = O(g2). 



(23) 



(24) 



The sum X^kk' IIkkKq ^ 0) is equal to p^K with p the 
average density and K the compressibility of the elec- 
tronic system. The latter provides a measure of the gap 
in the electronic spectrum because K vanishes as a func- 
tion of the chemical potential in the gap region. 

Equations (PT1) - ([M)) . respectively, can be considered as 
an orbital resolved closed form to describe the metal- 
insulator transition in terms of the polarizability or the 
compressibility of the electronic system, respectively. 
This is a primary tool to characterize the correspond- 
ing ground state. Such orbital selective sum rules are 
particularly useful in the proximity to a Mott insulating 
phase where different from the band insulator the internal 
degrees of freedom, orbital (and spin), approximatively 
survive. In order to fulfil the sum rules above, in case 
of the insulating state, in contrast to the metallic state, 
off-diagonal elements of the polarizability matrix describ- 
ing nonlocal polarization processes necessarily must oc- 
cur and interfere in order to correlate the charge response 
in such a way that Eqs. ([^5]) and ([M]) are satisfied. 

Physically this is related to the fact that in an insu- 
lator a perturbation, e.g., a change of the electron- ion 
potential, is only incompletely screened. Consequently, 
the selfconsistent change of the potential at the orbital 
degree of freedom n' nonlocally contributes to the CF's 
at the orbital k 7^ /t' in the unit cell. This is quite differ- 
ent from the metallic case where the diagonal elements 
of Hkk' dominate and off-diagonal elements can be ne- 
glected in most cases, i.e., only the PDOS Z^{£p) con- 
tribute. 

Approaching the delocalization-localization transition 
from the metallic region when p-doping is decreased ac- 
cording to our modeling of the cuprates the CuSd com- 
ponent of the wave function is admitted to become in- 
compressible, insulator-like in the underdoped state, i.e., 
in terms of the sum rule we have^ 

r 0{q) , Cu3d 
^H,,,(q^O)= ^ (25) 
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for the p-doped materials. In general, an orbital of type 
K is defined to be compressible or incompressible, respec- 
tively, in case the sum rule in Eq. ()2ip for a metal or 
Eq. (|23p for an insulator is satisfied for the orbital in 
question. 

In contrast to the p-doped case we propose in this work 
for the ri-dopcd materials in the underdoped state 

r 0{q) , 02p 

En«K'(q-o) = <^ ^ , (26) 

k' [ Zfi{ep) , otherwise 

i.e., the 02p component becomes incompressible, local- 
ized while the CuSd and Cu4s component remains metal- 
lic, delocalized as in the optimally doped metallic phase 
but with a renormalized density of states Zk(£f). We 
call a state defined by the sum rules from Eqs. (|25p and 
(j26p , respectively, a " strange metaF , because not all or- 
bitals in the system are compressible as in the normal 
metallic state related to the sum rule from Eqs. (PT|) . 
(|22p representing in general a Fermi liquid. 

Nevertheless, the total compressibility is never zero in 
such a strange metallic state and a real space organiza- 
tion for the low lying charge excitations is achieved in a 
complementary way in both cases. Quite generally, the 
incompressible regions compete with overall metallic be- 
haviour and with superconductivity. In our modeling of 
the underdoped phase of the HTSC's the PDOS for Cu3d 
is suppressed at the Fermi level in case of p-typc HTSC's 
and for 02p in the case of n-type materials. So, we have 
for p-doped materials an orbital selective compressible, 
metallic charge response for the holes in the 02p orbitals 
with a renormalized PDOS. For n-doped materials the 
charge response is metallic at the Cu sites but nonmetal- 
lic for the 02p orbitals simulating the absence of a hole 
contribution in the underdoped non-superconducting ma- 
terial. Altogether, we have a sudden loss in the density of 
states (pseudogap) at the Fermi level in both cases. Thus, 
the quasiparticle picture {Zfi^ep) 7^ for all k) consistent 
with our modeling of the metallic-, superconducting state 
from optimal to overdoping no longer holds for the under- 
doped regime. In the latter, we have a renormalized den- 
sity of states at ep which most likely is due to an orbital- 
dependent distortion of the bandstructure and a corre- 
sponding reconstruction of the large FS to some smaller 
one. Finally, in the optimally n-doped metallic state the 
formerly incompressible 02p orbital is allowed to become 
compressible metallic too as in p-doped LaCuO and holes 
can now contribute to the metallic properties and super- 
conductivity in n-typed materials consistent with the de- 
velopment of an LDA-like FS of hole type in NdCuO. In 
summary, in our approach the different ground states of 
the cuprates are distinguished by certain orbital resolved 
incompressibility- compressibility transitions related to 
corresponding transformations of the FS. 

Off-site CF's at the center of the planar oxygen squares 
can be introduced additionally, simulating specific cova- 
lence effects as described in Ref. [gj in order to improve 
selectively the phonon branch with the scissor mode at 



the X point. In this mode the in-plane oxygens vibrate 
perpendicular to the CuO bonds. These specific off- 
site CF's virtually have no effect on the other phonon 
branches. 

By our modeling of the underdoped state in both cases 
a partial ordering of the conducting carriers in the CuO 
plane is obtained. With regard to our modeling of the 
metallic state for optimally to ovcrdopcd HTSC's such 
an ergodicity breaking delocalization-localization transi- 
tion in the underdoped state in terms of the compressibil- 
ity means a compartmentalization of configuration space. 
This means that some parts of direct space cannot be ap- 
proached or are hardly accessible to the charge carriers 
in the pseudogap state. On the other hand, this should 
be accompanied by a reciprocal compartmentalization of 
momentum space, according to the uncertainity relation- 
ship Ax Ak ~ 1. Indeed, the ARPES studies in the 
underdoped cuprates, Refs. [25ll26ll27ll28l . point to such 
a compartmentalization of reciprocal space. The corre- 
sponding transformation of a large FS into Fermi surface 
arcs or small Fermi surface pockets, respectively, would 
result in a loss of density of states at the Fermi level 
(pseudogap) consistent with our modeling. 

III. PHONON DISPERSION OF NDCUO IN 
THE INSULATING AND METALLIC PHASE 

A. Phonon dispersion in the insulating state 

1. The ionic reference system - RIM 

For a definitive discussion of the phonon renormaliza- 
tion induced by the nonlocal EPI effects of CF and DF 
type, mediated by the second term in the dynamical ma- 
trix from Eq. ([5]), a quantitative reference model for 
the calculation of the phonon dispersion based essentially 
on the ionic component of binding is necessary. Such a 
model sketched in Sec. |TT1 representing approximately 
the local EPI effects, is provided by the ab initio rigid- 
ion model, extended via covalent ion softening and scal- 
ing of the short-ranged part of certain pair potentials. 
The results for the phonon dispersion of NdCuO along 
the main symmetry directions A (1,0, 0), E ^ (1,1,0) 
and A ~ (0, 0, 1) are shown in Fig. [T]where it is compared 
with the experimental results from inelastic neutron scat- 
tering (INS) for the insulating stated. The static effec- 
tive charges found for the model are Nd^'^'^+, Cu^'^^+, 
0^'^^~, Ol'^'^~ , which should be compared with the nom- 
inal ionic charges Nd^+, Cu^+, 0^~, Ol~ . The difference 
between these two sets of charges indicates the degree 
of ion softening. We see from Fig. [T] that the calcu- 
lated curves partially represent the experimental results, 
already quite well. However, the high-frequency opti- 
cal phonon modes (e.g. the Eu-modes at the F point) 
are considerably overestimated. This points to a miss- 
ing electronic polarization mechanism in the RIM. The 
strong ionic component of binding indicates that non- 
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FIG. 1: Calculated phonon dispersion of NdCuO in the 
main symmetry directions A ~ (1,0,0), E ^ (1,1,0), and 
A ~ (0, 0, 1) as obtained for the rigid-ion model explained in 
the text. The experimental data for insulating (undoped) 
NdCuO are taken from Ref. [39l . The diverse symbols 
representing the experimental results from inelastic neutron 
scattering (INS) mark the different irreducible representa- 
tions (ID). The arrangement of the panels from left to right 
is as follows: A1IA2 (• ■ ■ , o), A4 ( — , □)|A3|Si|E2 (• ■ ■ , o), 
S4 ( — , □)|E3|Ai (— , □), A2 (■ ■ • , o)|A3|. The frequencies 
are in units of THz. 



local DF's should be an important electronic degree of 
freedom in the charge response. 

We note that the RIM leading to the dispersion in Fig. 
[T] also gives reasonable structural data. So we obtain 
for the energy-minimized structure for the planar lattice 
constant a = 4. 155 A, for the lattice constant along the 
c-axis c = 12.741 A and for the internal position of the 
Nd ion .z(Nd) = 0.153 in units of c. The experimental 
values are a = 3.9451, c = 12.17lA and z(Nd) = 0.149, 
Ref. H 

From our numerical calculations for a purely ionic RIM 
with nominal ionic charges, not shown here, we find the 
following effects of covalent softening and scaling on the 
phonon dispersion. Scaling of the short-ranged part of 
the pair potentials has the qualitative effect to stabilize 
all formerly unstable modes found in a purely ionic RIM. 
Moreover, we find a decrease of the width of the spectrum 
towards the width of the experimental result. Ion soften- 
ing also decreases the spectral width but hardly stabilizes 
the unstable modes. Altogether, a considerable improve- 
ment of the calculated dispersion is achieved as compared 
to the results from a purely ionic model using nominal 
ionic changes and neglecting scaling. 



2. RIM plus dipole fluctuations 

Likewise as in our calculations for the p-doped materi- 
als LaCu02i, YBaCuO^, Bi-2201 and Bi-2213i^, a sig- 
nificant improvement of the calculated phonon dispersion 
as compared to the RIM can be obtained also for NdCuO 
by introducing additionally DF's into the modeling, see 





Cu 


Oxy 


0^ 


Nd 


ab-initio 


8.9 


7.2 


8.2 


12.5 


xy 


30% 


30% 


40% 


40% 
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40% 


100% 


40% 


40% 



TABLE I: The calculated result for the dipole polarizability 
according to the Sternheimer method^ in units of a% is given 
in the first row for Cu^-^2+, 0^^^", Oi'^*" and Nd'^-^^+ . The 
anisotropic reduction of the polarizability for the insulator is 
presented in the last two lines in percent. 100 percent stands 
for the ab initio result. 
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FIG. 2; Same as in Fig. [T]with the calculated results obtained 
taking additionally anisotropic dipole fluctuations (DF's) into 
account. 



Fig. [2] As in the calculations for the p-doped materials 
the dipole polarizability turns out to be anisotropic for 
the different ions, i.e., the calculated ab initio values of 
the dipole polarizability a as obtained from the Stern- 
heimer method for the isolated ions have in general to 
be reduced in a suitable way in the crystal. The actual 
values used in the calculations are listed in Table [J .This 
leads to a better agreement of the phonon dispersion with 
the experiments. A comparison of the calculated data 
in Fig. [2] with the experimental ones shows an overall 
good agreement. This means that already the allowance 
of DF's in the electronic screening process leads to rea- 
sonable results for the phonon dispersion of insulating 
NdCuO reflecting its ionic nature. 

Contrasting the calculated results of Fig. [T]with Fig. [2] 
the important effect of the DF's on the disperion can be 
figured out, in particular for the high-frequency optical 
modes. For example, the frequencies of the longitudinal 
E„ modes as well as the corresponding LO-TO splittings 
are properly decreased. The latter come out too large in 
the RIM, but now are in agreement with the experiments. 
As a minor deficit the width of the spectrum is still a 
bit too large. The two Ai-branches with the highest 
frequencies overestimate the frequencies at the F point 
(Ell modes) by about 0.5 THz and the highest Ei branch 
is about 1 THz too high in frequency. 
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FIG. 3: Same as in Fig. [T]with the calculated results achieved 
taking additionally DF's as in Fig. [2] and charge fluctuations 
(CF's) according to earlier calculations for LaCuO, Ref. [2^ . 
into account. 



FIG. 4: Same as in Fig. [T]with the calculated results for insu- 
lating NdCuO received taking additionally DF's as in Fig. [2] 
and optimized CF's with a strongly reduced CF-polarizability 
of the 02p orbital for NdCuO into account. 



3. RIM plus dipole- and charge fluctuations 

In a further step besides DF's as in Fig. [2] CF's are 
additionally introduced into the modehng according to 
the sum rule from Eq. for the charge response of an 
insulator. As in the case of LaCuO, Ref. [2J, CF's arc 
allowed at the Cu3c? and 02p orbitals in the CuO plane. 

In a first attempt we adopt the matrix elements of the 
electronic polarizability IIkk' foi' the CF's from our cal- 
culations for LaCuO^i, leading in the latter case to a 
good representation of the INS data. The correspond- 
ing results for NdCuO of such a modeling arc displayed 
in Fig. [3l As compared with Fig. [2j we recognize in 
contrast to the experiment a very strong softening of the 
dispersion of the two characteristic Ai and Si branches 
which display the anomalous behaviour widely seen in 
the metallic phase of the HTSC's, see below. In partic- 
ular, at the F-point a very large decrease in frequency 
is obtained for the high-frequency Eu-modcs leading to 
a strong disagreement with the experiment. From our 
numerical studies we find that the reason for this dis- 
agreement is the large matrix elements for the electronic 
polarizability Il{0xy2p) and 'n{0x2p — Oy2p) of the 02p 
orbital adopted from the case of LaCuO. 

In order to obtain the result for the dispersion of in- 
sulating NdCuO including DF's and CF's displayed in 
Fig. [31 which is in good agreement with the experiments, 
the polarizability of the 02p orbital concerning CF's has 
to be reduced to practically zero in NdCuO. Any en- 
hancement of the 02p-polarizability leads to a charac- 
teristic softening of the two Ai and Si branches with 
highest frequency towards the F-point in disagreement 
with experiment. The corresponding high-frequency E„ 
modes at the F-point are of CuO bond-stretching type, 
see Fig. [Sj and are consequently very sensible to the 
CF-polarizability of the 02p orbital. Thus, we find from 
our phonon calculations a significant different behaviour 
of the 02p orbital in the charge response in terms of 
CF's when comparing LaCuO with NdCuO. The nearly 
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FIG. 5: Displacement patterns 
of the two E„ modes of CuO 
bond-stretching type at the F 
point with the highest frequen- 
cies for the model from Fig. 
m The frequency of the modes 
given below the patterns is in 
units of THz. 



vanishing CF-polarizability of 02p in NdCuO points to 
stronger localization of that orbital in NdCuO than in 
case of LaCuO and is an expression of an enhancement 
of the ionic component of binding. 



B. Phonon dispersion in the metallic state 

In order to describe the kinetic part of the charge 
response for the metallic optimally doped (supercon- 
ducting) phase of NdCuO we apply for the modeling 
of the sum rule for the metallic state from Eq. 
(j2ip . Likewise as in our calculation for optimally p-doped 
LaCuO^i, CF's are only permitted in the CuO plane. In 
detail there are CF's of the CuSd and 02p orbitals ad- 
mitted as in the insulating state. Moreover, CF's of the 
delocalized Cu4s orbital are additionally allowed. The 
latter, introduced via the n(Cu4s) matrix element of the 
polarizability, have been found in the p-doped materials 
to be important for a correct description of the character- 
istic phonon anomalies in the highest Ai and Si branch, 
Rcfs. [Ullllllllllil. The growing importance of the 
Cu4s component in the wavcfunction, which leads to en- 
hanced delocalized CF's at the Cu sites in the screening 
process, has been related to a corresponding softening of 
the anomalous OBSM as compared to calculations where 
the Cu4s degree of freedom is neglected, see e.g. Ref. IH. 
So, from these calculations additionally to Cu3d and 02p 
an increasing contribution of the delocalized Cu4s com- 
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poncnt to the metallic state of the HTSC's is important 
to understand the charge response and the corresponding 
phonon softening of the OBSM, at least at higher doping 
levels. 

The importance of the Cu4s orbital for a realistic de- 
scription of the electronic structure of the HTSC's also 
has been pointed out in Refs. HIZsl . In this work, the 
hopping range has been identified as an essential material 
dependent parameter and the intralayer hopping beyond 
nearest neighbours as well as interlayer hopping proceeds 
via the Cu4s orbital. It is further concluded that materi- 
als with higher Tc.max have larger hopping ranges and in 
materials with highest Tc.max the axial orbital, essentially 
a hybrid between Cu4s, Cu3dz2_]^ and apical oxygen 2pz 
ist almost pure Cu4s. Moreover, the importance of Cu4s 
for an accurate description of partial charge distributions 
in the HTSC's is pointed out in Refs. 

In Fig. [5] we present our prediction of the phonon dis- 
persion for metallic NdCuO. We obtain a good agreement 
with the experimental results as measured so far by INS, 
Ref. [isl . In particular, the phonon anomalies are well 
represented. Our calculations support the INS results 
which on her part to some extent deviate from the cor- 
responding inelastic x-ray scattering (INX) results, Refs. 
[igIIitI . in particular as far as the phonon anomalies are 
concerned. The INS experiments find the anomalous oxy- 
gen half-breathing mode (Ai/2 anomaly) at ^(0.5,0,0) 
with a frequency of 12.5 THz, which is close to our calcu- 
lated value of 12.868 THz, while in the INX results this 
mode is preliminary identified at ^(0.2, 0, 0) with a fre- 
quency of about 12 THz. Our calculated data in Fig. [6] 
for the two Si branches with the highest frequencies arc 
also in very good agreement with the INS data while the 
INX results fall about 1.5 THz below the neutron data 
for the lower of the two Ei branches with the anomalous 
planar oxygen breathing mode (Og ) at the X point of 
the BZ. 

The calculated results shown in Fig. [6] for the metal- 
lic phase of NdCuO have been obtained with a model 
where the matrix elements H„„ have been increased for 
the Cu4s orbital but strongly decreased for 02p com- 
pared to our results for metallic LaCuO, Ref. [2l|. Thus, 
the tendency to a stronger localization correlated with 
a much smaller CF-polarizability of the 02p orbitals in 
NdCuO, i.e. a n-doped material, as compared to LaCuO, 
i.e. ap-doped material, as found already in the insulating 
state ist also present in the metallic state. 

An enhanced localization of 02p in NdCuO tends to 
stabilize the antiferromagnetic spin correlations in com- 
parison to a p-dopcd material like LaCuO, because the 
localization of the 02p orbital is favourable for superex- 
change. Such a trend is consistent with the fact that the 
phase diagram of the cuprates is asymmetric with respect 
to electron and hole doping. For the n-type materials the 
antiferromagnetic phase extends much further with dop- 
ing and competes with superconductivity. On the other 
hand, allowing the 02p orbitals to become compress- 
ible, metallic with a significantly larger CF-polarizability, 
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FIG. 6; Same as in Fig. [T] with the calculated results for 
metallic NdCuO taking additionally DF's and CF's as dis- 
cussed in the text into account. The experimental INS re- 
sults for optimally doped NdCuO shown by the diverse sym- 
bols, are taken from Ref. [l^ . The three black dots (•) in the 
figure indicate the anomalous high-frequency oxygen bond- 
stretching modes (OBSM) Ai/2, Of and Of , respectively. 



i.e. more delocalizcd as in our modeling of the p-doped 
HTSC's, makes the whole system more metallic, reduces 
antiferromagnetic fiuctuations and triggers superconduc- 
tivity already at a low doping level. 

The enhanced Cu4s CF-polarizability and the de- 
creased 02p CF-polarizability in a 7i-doped material like 
NdCuO as compared to a p-doped one seems to be con- 
sistent with an electron-hole asymmetry introduced by 
the doping process. While in the n-doped case materi- 
als electrons go to the hybridized orbitals of Cu3c? and 
Cu4s character during p-doping the holes are generated 
in the 02p orbitals and enhance the CF-polarizability of 
the latter as suggested by our calculations of the phonon 
dispersion. 

Altogether, our calculations of the phonon dynamics 
inp-doped and n-doped HTSC's point to a very different 
magnitude of the CF-polarizability of the oxygen ions in 
these materials. Thus, one might speculate that the CF- 
polarizability of oxygen should be indirectly relevant for 
both, antiferromagnetic behaviour and superconductiv- 
ity. The related stronger localization of the 02p orbital 
seems to favour in the n-doped material, via enhanced 
antiferromagnetic spin correlations, the competition be- 
tween antiferromagnetic order and superconductivity. 

As far as the modeling of the DF's is concerned the 
anisotropic dipole polarizability is modified when com- 
pared with the insulator. The dipole polarization in x- 
and j/-direction of the Cu and O^y ions in the CuO plane 
is assumed to be zero because of the metallic screening 
in the plane. Such a modeling already has been applied 
successfully for LaCuO in Ref. [l^. Moreover, we have 
reduced the dipole polarization of the apex oxygen in z- 
direction in the metallic state. This would be consistent 
with a better screening in c-direction via the delocalizcd 
Cu4s charge fluctuations, additionally allowed. The ac- 
tual values of the dipole polarizability arc given in Table 
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100% 


40% 


z 
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100% 


25% 


40% 



TABLE II: The calculated result for the dipole polarizability 
according to the Sternheimer method^ in units of a| is given 
in the first row for Cu^-^^+, O;^^^", Oj'^^' and Nd2-^^+ . 
The anisotropic reduction of the polarizability for the metal 
is presented in the last two lines in percent. The — denotes 
zero polarizability. 100 percent stands for the ab initio result. 

m 



C. PHONON DISPERSION ACROSS THE 
INSULATOR-METAL TRANSITION 

In this section wc propose a modeling of the insulator- 
metal transition in NdCuO via the underdoped phase 
characterized by the sum rule for the charge response 
given in Eq. (|26p . Looking from the perspective of the in- 
sulator this means that in such localization-delocalization 
transition expressed by the sum rule the 02p orbital re- 
mains incompressible, insulator- like (localized), while the 
CuSd and Cu4s orbital becomes compressible, metallic 
upon n-doping. 

In Fig. [7]we present a sequence of model calculations of 
the high-energy phonon modes related to the Ai, Si and 
Ai branches, respectively. The calculations include the 
anomalous OBSM i.e. the oxygen half-breathing mode 
(Ai/2 anomaly), the planar oxgen breathing mode at 
the X point (Og ) and the axial oxygen breathing mode 
Oj, at the Z point of the BZ, see Sec. IIVI for the cor- 
responding displacement patterns. In the leftmost panel 
the results for the insulator are shown and in the right- 
most panel those for the metal, compare with Fig. 2] and 
Fig. [6l respectively. 

For the description of the change of the phonon disper- 
sion across the insulator-metal transition in a first step we 
replace the dipole polarizations of the insulator (model 
II) by those of the metal. The result of this fictitious 
system is shown in the second panel from the left (model 
UO). The system is still an insulator, because the po- 
larizability related to the CF's has not been changed in 
the model, i.e. the Cu3d and 02p orbital is still incom- 
pressible according to the sum rule in Eq. ([^5)) . The 
next four panels display our results for the phonon dis- 
persion across the transition via an underdoped strange 
metallic state. The corresponding sequence of models is 
denoted U1-U4 in Fig. [7] and the sum rule from Eq. 
is used for the modeling, i.e. the 02p orbital is kept in- 
compressible. In order to satisfy this requirement some 
of the matrix elements of the polarizability must be cho- 
sen from the outside because the system of equations is 
under-determined. In particular we take the important 
n(Cu3(i) and n(Cu4s) component of the matrix as free 
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FIG. 7: Prediction of the anomalous phonon branches of Ai, 
El and Ai symmetry, respectively, across the insulator-metal 
transition. The different models are denoted from left to right: 
II, UO, Ul, U2, U3, U4, Ml. I, U, M denote the insulating, 
underdoped and metallic state, respectively. The unit of the 
frequency is THz. A detailed description of the calculated 
results in the different panels is provided in the text. 



parameters in order to calculate their possible effect on 
the phonon dispersion when they are varied across the 
transition. 

The definite variation of these matrix elements for 
modeling the sequence U1-U4 can be extracted from Ta- 
ble [Till The partial density of states ^^(eF) for the Cu3d 
and Cu4s orbital on the right hand side of Eq. ([26]) 
is taken to be equal to n(Cu3(i) and n(Cu4s), respec- 
tively. This is approximatively true, because the diagonal 
elements of 11^^' dominate for the compressible, metal- 





11, uo 


Ul 


U2 


U3 


U4 


Ml 


Cu3d 


0.05 


0.05 


0.20 


1.40 


2.800 


2.800 


Cu4s 




0.02 


0.02 


0.04 


0.075 


0.075 


02p 


0.01 


0.01 


0.01 


0.01 


0.010 


0.030 


Off 


0.02 


0.02 


0.02 


0.02 


0.020 


0.020 



TABLE III: Diagonal components n(Cu3d), n(Cu4s), 
n(02p) and n(Off) of the polarizability matrix 11^^/ used 
in the different models for the calculation of the phonon dis- 
persion of NdCuO across the insulator-metal transition via 
the underdoped state. II denotes the insulator, Ml the metal 
and U0-U4 the underdoped state simulating different doping 
levels. Note the different dipole polarizations in II and UO as 
discussed in the text. Units for n^K' are in eV~^ . 





11, UO 


Ul, U2, U3, U4 


Ml 


0^2p-Oy2p 


-0.00100 


-0.00100 




Cu3d-Off 


-0.00725 


-0.00725 


-0.00250 


Cu3d-02p 


-0.00525 


0.00100 




02p-0ff 


0.00225 


0.00225 


-0.00250 


Cu3d-Cu4s 




0.02500 




Cu4s-02p 




-0.00625 





TABLE IV: Values of the off-diagonal elements of the CF- 
polarizability in units of eV~^ for all the models discussed in 
the text. n(Cu4s — Off) is virtually zero and not given in the 
table. 



lie components. The matrix elements Il{Oxy), n(OfF), 
n(0:, - Oy), n(Cu3d - Off), n(0^y - off) are taken 
from the model II of the insulator, where the Cu4s or- 
bital is not involved. Finally, the remaining four un- 
known parameters, i.e. n(Cu4s — Off), n(Cu3(i — Cu4s), 
n(Cu4s — O^y) and n(Cu3(i— O^y) are determined from 
the sum rule in Eq. . Of course such a procedure is 
not unique, see in this context also our discussion of a 
similar modeling for LaCuO in Ref. [H. The complete 
list of the definite values for the matrix elements of the 
CF-polarizability for all the models used in our calcula- 
tions is summarized in Table Hill and HVl 

The modification of the charge response of the rele- 
vant orbitals CuSc? and Cu4s in terms of the polarizabil- 
ity matrix simulating doping in our approach is given in 
Table IIIIl We assume that the main effect introduced by 
ri-doping is to populate the hybridized orbitals of CuSd 
and Cu4s character. Thus, the partial density of states 
Zk{£f) for these orbitals should increase and correspond- 
ingly the components n(Cu3(i) and n(Cu4s) of the polar- 
izability. In Table [ni] we specify the values for n(Cu3(i) 
and n(Cu4.s) which have been used for the calculation 
of the phonon dispersion in the models U1-U4 as shown 
in Fig. [71 To act as a reference the parameters for the 
insulating and metallic state, II and Ml, are also given. 
Increasing n(Cu3d) and n(Cu4s) in the sequence U1-U4 
simulates an increase of n-doping in our modeling. This 
is carried out until the magnitude of these quantities are 
conform with model Ml for the metallic phase. In the 
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latter, the formerly incompressible 02p orbital is allowed 
to become compressible and now additionally may con- 
tribute to the metallic properties in NdCuO. At the end, 
we obtain the phonon dispersion for the model Ml shown 
in the rightmost panel where all orbital degrees of free- 
dom are compressible, metallic. A good agreement with 
the experimental dispersion curves including all their un- 
usual features is achieved. 

Concerning the change of the electronic properties 
across the insulator-metal transition via the strange 
metallic state, which should display non-Fermi-liquid 
properties according to our modeling^!, the main effect 
is related to an abrupt change of the electronic structure 
near the Fermi level and is globally expressed by an in- 
crease with doping of the partial density of states Zi^lep) 
for the compressible, metallic Cu3d and Cu4s orbital, re- 
spectively. The sudden loss in the DOS at the Fermi level 
(pseudogap) when passing from the model of the normal 
Fermi liquid with a large FS to a model for the strange 
metallic state with some reconstructed smaller FS by an 
orbital resolved compressibility-incompressibility transi- 
tion can be expected to show up in all physical proper- 
ties which scale with the DOS at the Fermi level, like 
thermodynamic properties, e.g. the electronic specific 
heat coefficient, transport properties as the resistivity or 
charge- and spin correlation functions. In context with 
our generic modeling of the strange metallic state in the 
cuprates and the related FS reconstruction, most likely 
conspiring with orbital-dependent band distortions, it is 
interesting to note that quite recently the existence of a 
small FS has been approved experimentally by the obser- 
vation of quantum oscillations of the Hall resistance in an 
underdoped cuprate, see Ref. and also Ref. [4^. Fi- 
nally, such a doping dependent FS reconstruction could 
be a means to explain, at least partly, the doping de- 
pendence of the unusual normal-state properties of the 
cuprates. 

While the pseudogap in the normal state of the 
cuprates is due to the incompressible, insulator-like states 
in our modeling, a distinct gap below Tc related to su- 
perconductivity can result from pairing of the doped 
charge carriers at a lower energy scale in the compress- 
ible, metallic states. Interestingly, according to the sum 
rules in Eqs. (QSl), (ESI) the condition for incompress- 
ibility and compressibility are not independent because 
of the off-diagonal elements of the polarizability matrix. 
Thus, it can be expected that below Tc the pseudogap 
and the superconducting gap are correlated and above 
Tc the small FS (arcs or pockets) and the pseudogap 
are related too. For example, a negative value of the 
(Cu3(i-0xy) matrix element necessary to make the CnSd 
state incompressible according to Eq. (PSI) . as required 
for underdoping, simultaneously will decrease the PDOS 
for the compressible 02p state. This correlates with a 
decrease of the small FS in the underdoped normal state 
and predicts a reduction of the superconducting gap if a 
BCS scenario is assumed. On the other hand, the matrix 
elements of the polarizability for Cu3d and Oxy decrease 
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when doping is reduced towards the insulating phase. So, 
the pseudogap as a loss in the density of states at the 
Fermi energy will increase with less doping. 

As already mentioned, when going from the under- 
doped state via the optimally to the overdoped state 
in our modeling all the orbitals become compressible, 
metallic. The pseudopgap effect related to the local- 
ized, incompressible Cu3d orbital in underdoped p-type 
cuprates vanishes and the compressibility coming along 
with a delocalization of all states corresponds to a re- 
covering of the reconstructed small FS to a large one 
consistent with ARPES experiments. The formerly in- 
compressible Cu3d orbital now dominates the DOS at Ef 
in the modeling. Additionally, a (virtual) Cu4s contri- 
bution develops at ef, which has shown its important 
fingerprints in the forming of the OBSM phonon anoma- 
lies upon doping. Altogether, a critical point which sep- 
arates two qualitatively different electronic ground states 
is reached in our modeling when the incompressible or- 
bitals become compressible accompanied by enhanced 
fermionic particle-number fluctuations, a transformation 
of the FS and the vanishing of the pseudogap effect. 
The orbital resolved compressibility of Cu3d in p-type 
cuprates and of 02p in n-type cuprates is a certain mea- 
sure of the local particle number fluctuations in the cor- 
responding orbitals. The incompressibility of these or- 
bitals in the underdoped state corresponds only approx- 
imatively to locally conserved particle numbers, because 
of the finite interaction strength in a realistic approach. 
In case, such quantities would be strictly conserved the 
latter are linked by general principles with local gauge 
symmetry of the theory. Using reduced model Hamilto- 
nians (like the infinite U Anderson or Hubbard model) 
the gauge symmetry approach has been extensively dis- 
cussed in the literature in context with the cuprates, see. 
e.g. Ref. HO. Physically, the partial incompressibility ex- 
pressed by ^Cu3d(£_F) = for p-doped and ^o2p(£f) = 
for n-doped materials in our modeling and thus the cor- 
responding reduced local particle density fluctuations are 
mainly related to localization by strong correlations and 
localization effects of ionic origin, respectively. This con- 
strains the low-energy dynamics of the electrons in the 
underdoped state, is responsible for the destruction of 
the large FS of the Fermi-liquid and is responsible for 
the non-Fermi liquid behaviour, correlating with a re- 
constructed FS. 

A superconducting state may become possible through 
an attractive interaction leading to pairing and destabi- 
lization of the FS as a result of an interplay of spin, charge 
and lattice degrees of freedomi^,. Spin-fluctuations, 
phonons and in the case of the optimally to overdoped 
state also coherent CF's along the c-axis in form of low 
lying phonon-plasmon modes provide a retarded contri- 
bution to the pairing interaction, while CF's and resid- 
ual anti-ferromagnetic spin-fluctuations in the Cu-0- 
plane contribute a non retarded, approximatively instan- 
taneous, part to the interaction at a higher energy scale. 

The large FS in the optimally to overdoped state also 



gives rise to a corresponding larger phase space for pair- 
ing as compared with the underdoped phase, in case some 
small FS due to the localization of the Cu3d or 02p in 
p- or n-doped cuprates, respectively, is representative for 
this state. According to the additional compressibility 
of the Cu3d orbital or 02p, respectively, particle number 
fluctuations arc enhanced as compared with the under- 
doped, partial incompressible state. Thus, phase coher- 
ence is established more easily according to the particle- 
number-phase uncertainity relationship AiVA^ > 1- A 
similar duality also exists locally, Ref. [sil. On the other 
hand, the reduction of particle-number-fluctuations in 
the underdoped state may clear the way for the formation 
of electron pairs at a higher energy scale well above the 
onset of phase coherence at Tc. Such a two-particle con- 
tribution of preformed pairs would also be accompanied 
by a reduction in the DOS which would add to the loss 
in the DOS from the orbital selective incompressibility in 
the single-particle channel. 

Furthermore, it should be remarked that within our 
model for the p-doped cuprates according to (PS)) the 
compressible, metallic 02p component of the electronic 
state is always connected to superconductivity, while the 
CnSd and Cu4s component additional comes into play for 
optimal to ovcrdoping. On the other hand after the sum 
rule from (j26p the situation is reversed for the n-doped 
case, i.e. superconductivity in the underdoped phase is 
orbital selective in the spirit of our modeling. 

The incompressible, localized Cu3d orbital in p-doped 
and the 02p orbital in n-doped cuprates also can be ex- 
pected to promote antiferromagnetic spin fluctuations. 

An ab initio calculation of the polarizability matrix 
n for the HTSC's, say beyond LDA, which should in- 
clude the strong electron correlations and also the strong 
nonlocal EPI effects in the self energy of the electrons, 
seems not possible in the near future. In particular this 
holds for the underdoped state of the cuprates. Thus, 
the PDOS, Zf^{eF), Z„(eF) and the matrix elements of 
n appearing in the various sum rules (PT|) - (P5)) can be 
understood as quantities renormalized by these interac- 
tions. Their deflnite values are presently unknown. A 
LDA result possibly may be a guide for the optimally to 
overdoped phase. In general, the pseudogap and the su- 
perconducting gap and the correlation in our approach is 
of course affected by these renormalization efltects. In any 
case, in our modeling the matrix elements of 11 and the 
PDOS can be treated as microscopically well defined pa- 
rameters of the theory constrained by the rigorous com- 
pressibility sum rules. 

Finally, the strong doping dependent nonlocal EPI ef- 
fects accompanied by electronic CF's as found for the 
OBSM and, in particular, also for the nonadiabatic 
phonons of polaronic character (mixed phonon-plasmon 
modes) propagating in a small region around the c-axis of 
the HTSC'si^ are certainly important to understand the 
unusual doping dependent isotope ejfect in the cuprates, 
because of their contribution to the electronic self energy. 

As already mentioned earlier, n(02p) in model Ml is 
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much smaller as for the modeling of the optimally doped 
metallic state of p-doped LaCuO because of a stronger lo- 
calization of 02p in NdCuO. In case of LaCuO we have 
obtained n(02p) = 0.2eV-^, Ref. [H, which has to be 
compared with U{02p) = 0.03eV-^ in case of NdCuO 
(Table IIIip . This small value points to a much weaker 
pseudogap effect in a n-type material as NdCuO as com- 
pared with the p-doped cuprates. 

As remarked in the introduction, in the ARPES ex- 
periments a large LDA-like FS of hole type evolves at 
optimal doping in metallic (superconducting) NdCuO, 
similar as in p-doped LaCuO, from the electron pockets 
seen in non-superconducting underdoped NdCuO. This 
suggests that holes might play a similar role in both 
types of superconductors. Moreover, as discussed in Sec. 
IIIIBI the CF-polarizability of the delocalized Cu4s or- 
bital becomes important in the metallic phase pointing 
to a substantial mixing of Cu4s with Cu3d and 02p in 
the wavefunction at the Fermi level. Figuratively, this 
may be interpreted by the emergence of holes in the 02p 
states which creates room for a delocalization of the Cu 
related states in particular via the Cu4s orbital. The ef- 
fect of the insulator- (strange)-metal transition has as a 
consequence a qualitative renormalization of the phonon 
modes as can be extracted by comparing the results from 
the second and third panel of Fig. \7\ These qualitative 
changes arc due to the closing of the LO-TO splittings 
of the Ell modes along the A and E direction at the F 
point and the vanishing of the A2u discontinuities along 
the A-direction if an adiabatic charge response along the 
c-axis is assumed as in our modeling. 

From our calculations of the electronic charge response 
and the resulting phonon dynamics the following quali- 
tative physical picture is consistent with our modeling of 
the electronic state in the cuprates, see also the detailed 
discussion for the p-doped case in Ref. [2l|. In the under- 
doped p-type material metallic hole carriers are related 
according to the sum rule from Eq. (j25p to the compress- 
ible 02p orbitals in the CuO plane. Cu-sites are avoided 
by the holes because of the strong on-site Coulomb inter- 
action, Usd, of the CuSd orbitals in the cuprates which 
are treated as incompressible, insulator like. Thus, there 
will be no repulsive core of the Coulomb interaction at the 
Cu-site for the holes in the superconducting state. This 
is compatible for example with d-wave superconductivity 
which on his part is favoured by a specific feature of the 
cuprates, namely, their common CuO square lattice. 

In the optimally p-doped state the Cu3d/As orbitals 
also become compressible, metallic in our modeling be- 
cause doping of more holes in the CuO plane creates room 
for a delocalization and a corresponding gain in kinetic 
energy for the Cu related states, in particular for Cu4s 
by hybrization with the 02p orbitals. The mixing in of 
the Cu4s state also may contribute a weak s-component 
for pairing. Moreover, the loss of density of states at 
ep (pseudogap) expressed by Eq. (^5)) . due to the in- 
compressible Cu states in the underdoped material, is 
restored upon doping and so these states additionally be- 



come available for pairing. 

In underdoped n-type materials Cu3d/4s orbitals are 
allowed to be compressible according to Eq. while 
02p remains incompressible. This means that we have 
an electron dominated transport in the normal state and 
a repulsive core of the Coulomb interaction at the Cu- 
site which is unfavourable for s-wave pairing of electrons 
because of the large on-site Coulomb interaction Usd, 
while holes in compressible, metallic 02p orbitals are not 
present unlike to the case in the underdoped state of the 
p-doped materials, see Eq. (|25p . However, s-wave pair- 
ing could become a possibility with increased doping be- 
cause of the small on-site Coulomb interaction IJ^s ■ Note 
in this context; that the importance of the Cu4s orbital 
for the charge response across the insulator-metal transi- 
tion via the underdoped state is growing in the modeling 
of the phonon dispersion in Fig. [T] I?-wave pairing via 
hole carriers is not allowed in the modeling for the un- 
derdoped state because the 02p orbital is incompressible. 
On the other hand, in the optimally doped case the lat- 
ter becomes compressible, metallic too and hole carriers 
in the 02p orbitals can emerge like in p-doped materi- 
als. Remember that the FS of optimall y dop ed NdCuO 
is of the hole-type and LDA-like, Refs. |25||26| . similar as 
for optimal p-doped LaCuO. In this way the repulsive 
core at Cu of the Coulomb interaction is avoided for the 
holes. Altogether, this leads to the conclusion that holes 
in the 02p states should play a similar role for d-wave 
pairing in both types of superconductors. A conclusion 
that holes are responsible for the superconductivity also 
in 71-doped cuprates quite recently has been drawn for the 
electron-doped superconductor Pr2_a;Ce2;Cu04 from the 
experimental side, i.e. by measurements of the resistivity 
and Hall angle as a function of doping and temperature, 
Ref. M- 

Nonadiabatic modifications to investigate the charge 
response around the " ionic" c-axis have been employed in 
Refs. [lUi Given the nearly two-dimensional electronic 
structure of the cuprates and thus a very weak interlayer 
coupling electron dynamics and phonon dynamics will 
be on the same time scale in a small region around the 
c-axis. This needs a nonadiabatic treatment with dynam- 
ical screening of the bare long-ranged Coulomb interac- 
tion and leads to a strong nonlocal, nonadiabatic EPI of 
polaronic character with phonon-plasmon mixing in the 
metallic state of the HTSC's. Such a strong coupling of 
the phonons to the electrons along the c-axis acts against 
a coherent interlayer hopping along this axis, helps on a 
confinement of the electrons in the CuO plane and may 
trigger a metallic to a nonmetallic crossover in the c-axis 
resistivity in case of a sufficiently low lying c-axis plas- 
mon. The energetic position of the plasmon is of course 
material specific and depends on the anisotropy and dop- 
ing level of the compound, Refs. [l9ll38l . 

Following the dispersion of the Ai modes in the metal- 
lic phase shown in the panels to the right of the second 
panel we find a continuous variation of the curves. In 
particular for the Ai branch with the second highest fre- 
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quencies we detect the development of a local minimum 
that can be identified with the Ai/2 anomaly found in 
the optimally doped metallic phase of NdCuO by INS. 
Also the Ai branch with the highest frequency develops 
a local minimum in full agreement with the experiments. 
The evolution to the two highest Si branches in Fig. [7| 
also converges to the experimental results for the opti- 
mally doped material, compare with Fig. (5] From these 
calculations it follows that a continuous enhancement of 
the polarizability related to CF's of the Cu3d and in par- 
ticular Cu4s orbital in the underdoped strange metallic 
state leads to the correct result for the optimally doped 
case. 



IV. PHONONANOMALIES AND 
ANTICROSSING BEHAVIOUR 



In the electron-doped material NdCuO the disper- 
sion of the branches with high frequency comprising the 
anomalous OBSM is more complex than in the hole- 
doped LaCuO. This becomes obvious when contrasting 
the results displayed in Fig. |4] and Fig. [61 wi th the anal- 
ogous results for LaCuO, see e.g. RefsnHIHIH- The 
reason is that in NdCuO there is another branch of the 
same Ai-symmetry which interacts with the branch con- 
taining the OBSM. Moreover, in the metallic phase (Fig. 
[S]) there is additionally a third branch of the same sym- 
metry in the A direction close-by in frequency. Thus, we 
will have anticrossing phenomena being absent in LaCuO 
because of the different crystal structure. Note, that a 
different lattice structure for LaCuO and NdCuO is also 
not in favour of an electron-hole symmetry. 

From the perspective of LaCuO (T-structure) the apex 
oxygen is shifted to a new position in NdCuO (T'- 
structure), sec Fig. [51 Here the bonding environment 
for the oxygen has changed as compared to LaCuO and 
vibrations with higher energies of the oxygen ions paral- 
lel to the CuO plane in the NdO "layer" become possible 
(Fig. ^ which may interact with the OBSM. In this 
figure also the calculated frequencies of the anomalous 
OBSM are listed for both, the metallic and the insulat- 
ing phase. From these data and from Fig. [Hand Fig. [Hlan 
anomalous softening of the OBSM across the insulator- 
metal transition can be extracted, similar to our findings 
for the p-dopcd cuprates. This, ultimately supports the 
generic nature of the phonon anomalies in p-doped and 
n-doped HTSC's. 

It is quite instructive to compare the calculated renor- 
malization of the OBSM due to strong nonlocal EPI of 
CF-type across the insulator-metal transition in NdCuO 
and LaCuO, respectively. Quantitatively we obtain for 
of, Ai/2 and Of the following softening in THz (Nd- 
CuO first number, LaCuO second number). O^: 2.233, 
2.539; Ai/2: 1.578,1.452; Of: 2.275,2.962. While the 
in plane polarized OBSM show a similar renormaliza- 
tion, the softening for is stronger and clearly larger 
in LaCuO. As far as the anomalous large lincwidth of 
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FIG. 8: Elementary cell 
of tetragonal LaCuO 
(T-structure, left) and 
NdCuO (T'-structure, 
right). •: Cu, o: O, •: 
La, Nd 



the of mode found in the experiments Refs. IsHolflsl is 
concerned we refer to our discussion of this mode in a 
nonadiabatic phonon-plasmon scenario ^^i'^^ . 

Importantly, strongly coupling modes, like Og and 
of, which involve momentum transfer between anti- 
nodal regions of the Fermi-surface should contribute via 
the strong nonlocal EPI to electron self-energy correc- 
tions in the (tt, 0, 0)-region and in particular to the anti- 
nodal pseudogap in the normal state. Moreover, there 
is an enhanced contribution to the self-energy in the su- 
perconducting state because of the density of state en- 
hancement in these regions due to the opening of a gap 
with d-wave symmetry, Ref. Issl While in a d-wave su- 
perconductor Og connects regions with different sign of 
the order parameter, O^ connects regions with the same 
sign. The Ai/2 anomaly on the other hand involves mo- 
mentum transfer between nodal regions and could be im- 
portant to understand the corresponding self-energy cor- 
rections, in particular the nodal kink observed in ARPES 
experiments, Refs. [llH. 

Inspection of the displacement patterns in Fig. [9l 
demonstrates that the anomalous OBSM Og , Ai/2 de- 
fined by oxygen breathing vibrations in the CuO plane 
always exhibit the lower frequency of the two modes with 
the same symmetry. The modes with the higher frequen- 
cies are essentially confined in the NdO "layer". In the 
axial oxygen breathing mode Of at the Z point, corre- 
sponding to the apex oxygen breathing mode in LaCuO, 
all the out of plane oxygens vibrate in phase against the 
CuO plane. Therefore, this mode may induce in the 
metallic state strong CF's between the CuO planes (in- 
terplane charge transfer), similar as in the case of LaCuO, 
see e.g. Ref. [2l| and Sec. [Vl 

Now, let us study the anticrossing effects of the high- 
energy modes in NdCuO incorporating the anomalous 
OBSM in the A and S direction. We restrict our discus- 
sion to the metallic state. In Fig. [TOlthe two interacting 
Si branches are displayed with a higher resolution than 
in Fig. [51 From Fig. [TUl the anticrossing point can be 
localized at about q ~ 0.4. This may be approved by 
the calculation of the displacement patterns of the corre- 
sponding eigenvectors. Typical for anticrossing the latter 
have changed their character (Fig. Ilip when going from 
q = 0.3 to q = 0.5. 

The situation for the three interacting high-frequency 
branches of Ai symmetry is more complex. These 
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14.162 17.061 12.868 15.282 15.023 

16.395 18.081 14.446 15.751 17.298 

FIG. 9: Displacement patterns of the anomalous OBSM Ob (left), Ai/2 (middle), Of (right) for NdCuO. The calculated 
frequencies in units of THz for the metallic phase are listed in the first row and those for the insulating phase in the second 
row below the patterns. 



18k - 

17 ; - 

15- 

15- 

X 14 - y 

/ - 

I 12 : / \ 

IQ - 
9 ^ - 
h F .... I .... I .... I .... I .... I .... I .... I .... I .... I .... - 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

q 

FIG. 10: Phonon dispersion of the two interacting Ei ~ 
^q(l, 1,0) branches from Fig. |6] displayed with a higher res- 
olution in order to study the anticrossing effect. 

branches are shown with a higher resolution in Fig. 1121 
From our calculations we find an anticrossing of the two 
highest Ai branches between q = 0.2 and q = 0.3. In the 
sector between q = 0.7 and q = 1.0 {Z point) we find a 
complex anticrossing scenario between three interacting 
Ai branches. The highest mode at g = 0.7 ends as the 
lowest mode at the Z point. The lowest mode at g = 0.7 
becomes the Of mode with the second highest frequency 
at the Z point, compare with Fig. 1131 where the displace- 
ment patterns of the relevant modes are given for the F 
point (q = 0), A/2 {q = 0.5) and the Z point {q = 1.0), 
respectively. Finally, the mode with the second highest 
frequency at q = 0.7 winds up to the highest frequency 
at Z. 



V. PHONON-INDUCED CHARGE RESPONSE 
OF THE OBSM 

In Fig. [T3]the result of the calculations of the phonon- 
induced redistribution of the charge density Sp is given 
for the anomalous OBSM in the metallic phase of Nd- 
CuO. Part (a) displays for the Ai/2 anomaly the total 
redistribution including, both, the local, rigid contribu- 



tion as obtained from the local EPI defined by the RIM 
in Sec. IIIII as well as the nonlocal, nonrigid part of the 
charge response related to the nonlocal EPI effects me- 
diated by the CF's. Part (b) gives exclusively nonlocal 
contributions of dp, i.e. 

<5p„(r, qa) = ^ <5C(qa)p« (r - R^), (27) 

with the CF's according to Eq. (|13p and the form factors 
for the CF's as in Eq. ([T]). The total phonon-induced 
charge redistribution can be obtained by adding to Eq. 
(P7|) the rigidly displayed unperturbed densities, i.e. 

Sp,{r,(ia) = 

^ {p° (r - [RS(qa) + <(qa)]) - p° (r - R«)} . 

(28) 

p'^ is the density of the unperturbed ion from Eq. ([1]) 
and the displacement of an ion in Eq. In the 

figures shown, Sp > (full lines) means that electrons 
are accumulated in the associated region of space. The 
broken lines on the other hand indicate that electrons 
are depleted. The rigid part of the charge response can 
clearly be seen in Fig. [Tib near the displaced ions, while 
the nonlocal, nonrigid part in Fig. 114b demonstrates that 
the moving oxygen O^; ions in the half breathing mode 
generate via nonlocal EPI changes of the potentials at the 
silent Cu and Oy ions resulting in corresponding CF's in 
form of a charge transfer within and between the CuOy 
chains. Note, that there are no changes of the transfer 
integral between d- and p-orbitals for the silent Cu and 
Oj, ions, nevertheless there is a charge transfer, nonlo- 
cally induced. In case of the Og mode, Fig. [Tib , the 
moving 0^;^ ions induce CF's at the silent Cu ion and 
we obtain an electronic charge transfer from that Cu ion 
where the CuO bonds are compressed to the Cu ion where 
the bonds are stretched. According to these calculations 
we have a strong coupling of the phonon modes to the 
charge dynamics. The nonlocal EPI leads to a (dynamic) 
charge ordering by CF's in the form of localized stripes of 
alternating sign in the CuO plane which are interacting 
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FIG. 11: Displacement pattern of the two interacting Ei-modes from Fig. [10] at q = 0.3, q = 0.4 and q — 0.5 in the region 
around the anticrossing point locahzed at about q = 0.4. The frequencies of the corresponding modes given below the patterns 
are in units of THz. 
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FIG. 12: Phonon dispersion of the three interacting Ai ~ 
(1,0,0) branches from Fig. [6] shown with a higher resolution 
in order to study normal anticrossing in the sector between 
q — 0.2 and q = 0.3 and complex anticrossing in the sector 
between q = 0.7 and q = 1.0. q = 0.0 represents the F point 
and q = 1.0 means the Z — ^(0, 0, 1) point in the BZ. 



with the lattice vibrations and with each other. In case 
of the Ai/2 anomaly (half breathing mode) the charge 
stripes arc directed along the x- or y-axis, respectively, 
and for Og (planar breathing mode) along the diagonals 
of the CuO plane. The charge patterns appear instan- 
taneously because of the adiabatic approximation used 
in the calculation which is sufficient for these modes. It 
should be remarked that qualitatively the same density 
redistributions are calculated for Ai/2 and Og in the 
metaUic phase of LaCuO, see Refs. HHHIlllii 

The nonlocal induced CF's for the Of mode in NdCuO 
is shown in Figs. [Tib , d. Here we have the situation that 
the displacement of the Oz ions (Fig. [9] and [T4H) in the 
ionic layers excites nonlocally changes of the potential 
felt by the electrons in the CuO plane which on her part 
are responsible for superconductivity. These facts visu- 
alize the importance of such c-axis polarized phonons for 
superconductivity via the electron-phonon mechanism. 

The nonlocal coupling effects are an expression of the 
strong component of the ionic binding along the c-axis 



in the HTSC's, and similar charge redistributions as in 
Figs. [T4b . d also have been calculated for the apex oxy- 
gen breathing mode in LaCuO^i. These long-ranged 
Coulomb coupling effects are very special to the HTSC's 
and would not be possible in a conventional metal or su- 
perconductor because of local-screening by a high-density 
electron gas. The strong softening of about 2.3 THz for 
the of mode when going from the insulating- to the 
metallic phase can physically be understood by compar- 
ing the phonon-induced charge rearrangements for Of as 
calculated in the metallic- (Figs. [Tib , d) and insulating 
state (Fig. ^Ejp). As shown in Fig. [5] and Fig. [TiB the 
Oz ions move in phase against the CuO layers, similar 
as the apex oxygen mode in LaCuO. Hence, because of 
the weak screening one can expect this vibration to in- 
duce CF's in the CuO planes as seen. However, these 
CF's are qualitatively different in the metallic- and insu- 
lating state, respectively, because of the energy gap for 
charge excitations in the insulator. This helps to under- 
stand the anomalous softening of the mode across the 
insulator-metal transition. As has been shown in Ref. [13 
the CF's SCfi{q(T) from Eq. are constrained in the in- 
sulator with CF's allowed at the Cu and Oxy sublattices 
according to the following sum rule: 

^<5Ck(A^) = 0. (29) 

K denotes the CF's in the CuO layer. A (0, 0, 1). Thus, 
Eq. ([^51) particularly holds at the Z point. In contrast to 
the constrained expressed by Eq. (j29p for the insulating 
state, which means that local charge neutrality of the 
cell is maintained under a perturbation due to Of , no 
such a restriction is present for the metallic state. Con- 
sequently, in the insulating state only intralayer charge 
transfer according to Eq. ([29)1 are allowed, see Fig. [TSb . 
On the other hand, in the metallic state Of induces 
CF's at Cu and Oxy of the same sign in the whole CuO 
layer (Fig. 114b . d). This finally makes an interlay er 
charge transfer possible, which appears instantaneously 
in the adiabatic approximation and provides an effec- 
tive screening mechanism for the long-range Coulomb- 
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FIG. 14: Contour plots of the phonon- 
induced charge density redistribution 
6p, according to Eqs. ((27l), ([28} for 
the OBSM in metalhc NdCuO. The 
units are 10~''/a|. FuU lines ( — ) 
mean that electrons are accumulated 
in the corresponding region of space 
and broken lines (- ■ ■ ■ ) represent re- 
gions where the electrons are pushed 
away. The arrows denote the dis- 
placement of the oxygen ions (a) total 
charge redistribution, i.e. Spr + Spn for 
the Ai/2 anomaly (b) nonlocal part 
5pn of the charge redistribution for the 
Ai/2 anomaly, (c), (d) Spn for Of . 
(e) 5p„ for Of. 



interaction and a corresponding softening of Of . As al- 
ready remarked the weak interlayer coupling in the very 
anisotropic cuprates most likely makes necessary a nona- 
diabatic calculation for modes in a small q-spaee region 
around the c-axis including O^. In this case we find 
in our calculations for LaCuO a nonadiabatic, insulator- 
like charge response crossing over to a coherent adiabatic 
metallic regime outside this regio n^^-^^ . Within such a 
treatment the massive line broadening of Of found ex- 
perimentally in LaCuO can be understood. According to 
Ref. [l^ the Of mode in metallic NdCuO is also not well 
defined indicating a large intrinsic linewidth which may 
be understood along the same lines as in LaCuO. Even- 



tually, the instantaneous interlayer charge transfer in the 
adiabatic approximation is replaced by a c-axis plasmon 
at the Z point in the nonadiabatic regime. 

Comparing the results in Fig. [T5b for the charge re- 
distribution due to the half-breathing mode for insulat- 
ing NdCuO with corresponding results of LaCuO we find 
that the stripe patterns directed in x- or y-direction, re- 
spectively, found in LaCuO is destroyed in NdCuO most 
likely due to the small CF-polarizability of the 02p or- 
bital in NdCuO. On the other hand, we find for Og qual- 
itatively the same diagonal stripe pattern as in LaCuO. 
On the whole, the charge redistribution of the OBSM in 
the insulating state are more localized around the ions as 
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FIG. 15: Same as in Fig. [T4]but for insulating NdCuO. (a) 5p^ for the Ai/2 anomaly, (b) 5pn for Of . (c) i5pn for Ob . 



in the metallic state. 



VI. PHONON DENSITY OF STATES 

In the last topic of this paper we present calculations 
for the phonon-density of states of NdCuO in the insulat- 
ing and metallic state according to our modeling. Com- 
paring the spectrum of the insulating state in Fig. 116b . 
with that of the metallic state in Fig. [TBb and the cal- 
culated results for the corresponding atom-resolved par- 
tial density of states (PDOS) in Fig. [T7] characteristic 
changes across the insulator-metal transition can be de- 
tected. From the calculated results we find across the 
transition a decrease of the width of the spectrum to- 
gether with a redistribution of spectral weight in partic- 
ular in the high-frequency part. The latter redistribu- 
tion is essentially due to the softening of the anomalous 
OBSM, denoted as A(Of ), B(Of ) and C(Ai/2) in Fig. 

Fig. [T7k demonstrates that the low frequency part up 
to about 6 THz of the spectrum is dominated by the vi- 
brations of the Nd ions. From the results for the PDOS 
of the Cu ions in Fig. 117b . which are important up to 
about 8 THz, we conclude that the effect of the transi- 
tion is marginal. However, Figs. 17c and 17d display the 
essential effect of the transition on the phonon spectrum 
in the high-frequency range. Similar as in the case of 
LaCu02i there is a characteristic shift to lower frequen- 
cies of the spectral weight for the high-frequency oxygen 
phonon modes across the transition, indicating again the 
strong clcctron-phonon interaction in the metallic state. 
Moreover, we do not have just a rigid shift but a char- 
acteristic rearrangement of the high-frequency part of 
the spectrum related to the development of the OBSM 
phonon anomalies in the metallic phase. 

The calculated phonon density of states can be com- 
pared with the inelastic-ncutron-scattcring results as ob- 
tained in Rcf. for metallic NdCuO. There are obvi- 



ous peaks at 3.1, 12.3, and 15.7 THz in the experiments 
which correlate very well with significant peaks in the 
calculated spectrum, shown in Fig. [TBb . From Fig. [T7| 
we extract, that the peak at 3.1 THz is mainly due to 
the Nd vibrations with only a small contribution from 
the other ions. The peak at 12.3 THz is dominated by 
the Oxy vibrations with a considerable contribution from 
O2. Finally, in case of the peak at 15.7 THz. both, Oxy 
and Oz participate with about equal strength. 

VII. SUMMARY AND CONCLUSIONS 

We have calculated for NdCuO complete phonon dis- 
persion curves, total- and atom-resolved density of states 
and local as well as nonlocal phonon-induccd charge re- 
distributions across the insulator-metal transition. We 
have used our microscopic modeling for the undoped 
insulating-, the underdoped strange metallic- and the 
optimally doped metallic state in terms of orbital selec- 
tive incompressibility-compressibility transitions, differ- 
ent from the p-dopcd case studied earlier. In particular, 
we have investigated the possibly generic phonon anoma- 
lies (OBSM) in the HTSC's and their complex anticross- 
ing behaviour in NdCuO. A good agreement with the 
corresponding experimental phonon dispersion measured 
by INS has been found. Our calculations support the 
assignments of the anomalies as provided by INS as com- 
pared to the INX results. 

In order to obtain the calculated results, both, long- 
ranged Coulomb interactions of ionic origin as well as 
short-ranged repulsive interactions of electronic origin to- 
gether with a sufficiently large set of orbital degrees of 
freedom are necessary. We conclude that a purely ionic 
model, similar as in case of LaCuO, leads to a consider- 
able overestimation of the width of the phonon spectrum 
and some unstable branches. Including covalcnt correc- 
tions a suitable reference model for the insulating phase 
of the HTSC's can be constructed that cures the disad- 
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FIG. 16: Phonon density of states for NdCuO as calculated for 
the insulating state (a) and the metallic state (b) , respectively. 
The characters A, B, C denote the OBSM Of , Of , and the 



Ai/2 anomaly, respectively. 



vantages of the purely ionic description. Starting from 
such a rigid reference model it has been shown that non- 
local, nonrigid electronic polarization processes via local- 
ized CF's and DF's are crucial for a correct specification 
of the density response and phonon dynamics. 

From our calculations some insight into the strange 
metallic state of the cuprates and an electron-hole asym- 
metry introduced by doping can be gained. We have 
found a characteristically different behaviour of the 02p 
orbital in the charge response when comparing p-doped 
LaCuO to n-doped NdCuO, respectively. In both, 
the insulating and metallic phase, the electronic CF- 
polarizability of 02p is strongly supressed in NdCuO as 
compared with LaCuO, indicating a stronger localization 
of the 02p orbitals in the n-dopcd material. This points 
to an enhancement of the ionic component of binding. 
Simultaneously, the stronger localization of 02p tends to 
stabilize the antiferromagnetic spin correlations in com- 
parison to a p-doped material like LaCuO, because a 
localization of the 02p orbital can be expexted to be 
favorable for superexchange. Note, that for the n-type 
cuprates the antiferromagnetic phase extends much fur- 
ther with doping and from our calculations one might 



speculate that the polarizability of oxygen should play an 
indirect role for antiferromagnetism and superconductiv- 
ity. 

In our modeling we found for NdCuO an increased ma- 
trix element for the Cu4s CF-polarizability in the metal- 
lic state as compared to LaCuO. An increased occupa- 
tion of this delocalized orbital in the hybridized metal- 
lic state at the Fermi energy, likewise as in LaCuO is 
very important to comprehend the softening behaviour 
of the OBSM, 0^-, Of-, Ai/2-anomaly, found across 
the insulator-metal transition via the underdopcd state. 
The enhanced Cu4s-polarizability in a n-doped material 
seems to be consistent with an electron-hole asymmetry 
introduced by the doping process. 

We have presented a qualitative physical picture of the 
electronic state in the cuprates consistent with our mod- 
eling of the charge response which on his part leads to 
a good agreement of the calculated phonon dispersion 
with the INS results, in particular for the generic phonon 
anomalies. These modes along with c-axis modes like Of 
from the nonadiabatic region most probably should con- 
tribute to electron self-energy corrections like kinks and 
pseudogaps seen in ARPES experiments as well as to 
the unusual transport properties along the c-axis in the 
cuprates. The definite calculations have shown that the 
full set of orbital degrees of freedom, i.e. Cu3(i/4s and 
02p is essential. Qualitatively s- and d-wave supercon- 
ductivity has been discussed within our modeling and it 
is conjectured that holes in the 02p states should play a 
similar role for d-wave superconductivity in both, p-type 
as well as n-type cuprate superconductors. 

From our calculation of the complex anticrossing sce- 
nario in NdCuO seen in the experiments but being absent 
in LaCuO we conclude that the anticrossing is related to 
the different crystal structure of these materials. 

Our calculations of the phonon-induced charge redistri- 
butions performed for the anomalous OBSM, Og , Ai/2, 
indicate that the nonlocal EPI of ionic origin leads to (dy- 
namic) charge ordering by CF's in the form of localized 
stripes of alternating sign in the CuO plane, minimiz- 
ing the energy of the system and decreasing the phonon 
frequencies considerably. 

The strong frequency renormalization of Of across the 
insulator-metal transition is attributed by our calcula- 
tion to an interlayer charge transfer in the metallic state 
not possible in the insulator where only intralayer charge 
transfer is allowed. 

Finally, there is, similar as in LaCuO, a characteristic 
shift to lower frequencies of the spectral weight for the 
high frequency oxygen vibrations in the metallic phase 
accompanied by a strong rearrangement of the high- 
frequency part of the phonon density of states. This is 
essentially related to the development of the anomalous 
OBSM in the metallic phase, due to a strong electron- 
phonon interaction, also in the n-doped compounds. 
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FIG. 17: Comparison of the atom resolved phonon density of states of NdCuO between the insulating state (• • • ) and the 
metallic state ( — ). 
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